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Abstract 

We take the first steps towards the complete computation of the thermalization 
rate of the supersymmetric particles involved in electroweak baryogenesis by 
computing the thermalization rate of the right-handed stop from the imainary 
part of the two-point Green function. We use improved propagators includ- 
ing resummation of hard thermal loops. The thermalization rate is computed 
at the one-loop level in the high temperature approximation as a function of 
Mi R (T). We also give an estimate for the magnitude of the two-loop contribu- 
tions which dominate the rate for small M^ R (T). If the stop is non-relativistic 
with Mf R (T) 3> T, thermalization takes place by decay and is very fast. 
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It is now commonly accepted that the generation of the baryon asymmetry dur- 
ing the electroweak phase transition Q requires some new physics at the weak scale. 
Threfore electroweak baryogenesis in the framework of the Minimal Supersymmetric 
Standard Model (MSSM) has attracted much attention in the past years, with par- 
ticular emphasis on the strength of the phase transition and the mechanism of 
baryon number generation |3], |], |5|] . It has recently been shown both analytically |7], Q 
and by lattice simulations || that the phase transition can be sufficiently strongly 
first order if the lightest stop is not much heavier than the top quark, the ratio of the 
vacuum expectation values of the two neutral Higgses tan (3 is smaller than ~ 4 and 
the lightest Higgs is lighter than about 85 GeV. 

Moreover, the MSSM may contain, besides the CKM matrix phase, new CP- 
violating phases in the soft supersymmetry breaking parameters associated with the 
stop mixing angle and with the gaugino and neutralino mass matrices. In the MSSM 
large values of the stop mixing angle are restricted in order to preserve a sufficiently 
strong first order electroweak phase transition. Therefore, an acceptable baryon asym- 
metry may only be generated through a delicate balance between the values of the 
different soft supersymmetry breaking parameters contributing to the stop mixing pa- 
rameter, and their associated CP- violating phases ||, |J. In the MSSM it turns out that 
the main contribution to the baryon asymmetry comes from charginos and neutralinos 
and the phase of the parameter /i must be larger than about 0.1 to be responsible for 
the generation of the observed baryon asymmetry [|], |^, |5| . If the strength of the elec- 
troweak phase transition is enhanced by the presence of some new degrees of freedom 
beyond the ones contained in the MSSM, e.g. some extra standard model gauge sin- 
glets, light stops (predominantly the right-handed ones) and charginos/neutralinos are 
expected to give quantitatively the same contribution to the final baryon asymmetry. 

Extra and sizeable CP-violating phases in the MSSM are therefore a necessary 
ingredient for a successful electroweak baryogenesis scenario. If the particles involved 
in the process of baryon number generation thermalize rapidly, CP-violating sources 
however loose their coherence and are diminished. This phenomenon can be intuitively 
understood by means of the following example: let us focus on the right-handed stop 
current and imagine that a right-handed stop scatters off the advancing bubble wall and 
is transformed into a left-handed stop. If the latter scatters off the Higgs background 
once again, a right-handed stop reappears in the plasma; since in both interactions 
with the wall CP is violated at the vertices because of the explicit phase in A t fi, 
quantum interference may give rise to a final right-handed current. This, however, 
only takes place if, along their way, the stops do not interact with the surrounding 
plasma and disappear: large thermalization rates reduce the final baryon asymmetry. 

Moreover, CP-violating currents in supersymmetric baryogenesis are more easily 
built up if the degrees of freedom in the stop and in the gaugino/neutralino sectors 



are nearly degenerate in mass [|J. For instance, phases of fi smaller than 0.1 are only 
consistent with the observed baryon asymmetry for values of \fi\ of the order of the 
gaugino mass parameters. This is due to a large enhancement of the computed baryon 
asymmetry for these values of the parameters. This resonant behaviour is associated 
with the possibility of absorption (or emission) of Higgs quanta by the propagating 
super symmetric particles. For momenta of the order of the critical temperature, this 
can only take place when, for instance, the Higgsinos and gauginos do not differ too 
much in mass. By using the Uncertainty Principle, it is easy to understand that the 
width of this resonance is expected to be proportional to the thermalization rate of 
the particles giving rise to the baryon asymmetry [Q. 

It is therefore clear that the computation of the thermalization rate of the particles 
responsible for supersymmetric baryogenesis represents a necessary step towards the 
final computation of the baryon number. Despite its relevance, no substantial effort 
has been devoted to a detailed computation of the decay width of SUSY particles in 
the thermal bath. The goal of this paper is to take the first step towards a complete 
evaluation of the thermalization rate of the supersymmetric particles involved in the 
generation of the baryon asymmetry, right-handed stops, charginos and neutralinos. 
At present, we restrict ourselves to the computation of the thermalization rate Tr of 
light right-handed stops from the imaginary part of the two-point Green function in the 
unbroken phase of the MSSM. We use improved propagators including resummation 
of hard thermal loop; the thermalization rate is computed at the one-loop level in the 
high temperature approximation and an estimate is given for two-loop contributions. 
Finite mass effects as well as the computation of the thermalization rate for charginos 
and neutralinos will be presented in a longer publication |]TT |. 



Let us first present some technical details needed for the computation of the right- 
handed stop thermalization rate. A great deal of the formalism may be found in 



Cqj , but we briefly summarise here some details and also present some new necessary 
technical tools. 

The thermalization rate T# depends upon the imaginary part of the two point 
Green function E# via the relation 



(1) 



where uj = uj(k) is the energy of a given mode. In the present paper we shall focus on 
the long wave-length stops with uj ~ Mf . At the one- loop level the quantity Im E# 
receives different contributions from diagrams involving SM fermions, charginos, neu- 
tralinos, scalar and gauge bosons^. At two loop level there are several potentially 
relevant sources that contribute to the thermalization rate. The strong interactions 

x We assume R-parity conservation so that one of the two particles running in the loop is always 
a superpartner. 



contribute, due to diagrams including gluon interactions being proportional to the 
square of the strong coupling constant. Also, there are possible interactions of two 
right-handed stops appearing from D-terms of the superpotential. They are, how- 
ever, proportional to the fourth power of the electroweak gauge couplings (and so 
are two-loop diagrams containing electroweak gauge bosons) and therefore can be ne- 
glected. Possible interactions arise also from the F-terms of the superpotential. Of 
particular importance is the one proportional to the top-quark Yukawa coupling given 
by ht\H 2 \ 2 \tii\ 2 . The diagram including (only) this interaction is not suppressed by 
the coupling, but is expected to be small due to phase space suppression unless the 
right-handed stop is light. The very same argument can be used in the case of gluon 
interactions, too. 

The imaginary parts are cuts across the relevant diagrams and correspond to the 
differences between the absorption and the emission rates, or between decay and inverse 
decay rates. For brevity we denote these respectively by "absorption" and "decay". 
For all the fermions and bosons unbroken phase thermal corrections must be taken into 
account. They merely change the pole structure of the propagators and, to leading 
order in the temperature T and for the massless SM fermions, they are the hard 
thermal loops. Accordingly, there appear energy thresholds which, in particular for the 
fermionic loops, are rather complicated because of the complicated nature of fermionic 



dispersion relations, which have two branches called particles and holes |13| . However, 
the vertex corrections relevant in the present discussion contains no hard thermal loops 



and can therefore be omitted 14 



The plasma corrected left-handed Dirac fermion propagator reads generically (P M 
(p , p), andp= |p|) 



S(P) - P L p 2 ^J^PR, 



(2) 



where Pl and Pr are the left- and right-handed projections and fi is the bare mass, 
and 

7m F" = [1 + a{P)]P,Y + 6(P) 7 °. (3) 



In the case of a left-handed field, at high temperature T 3> /i the functions a and b 
are given by [|13 
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where is the plasma mass of the left-handed particle. An identical formula can be 
written for right-handed fermion component with replaced by m^. 



Eqs. (|4],|5|) apply to all particles with masses <C T, that is, to all the SM fermions 
in the unbroken phase. For Dirac fermions with non-zero bare mass m, the situation is 
more complicated because the left- and right-handed components couple. The Majo- 
rana propagators for massive neutralinos would be even more complicated [13]. In the 
present paper we shall focus on the limit T ^> m, which we believe cover most of the 
cases of interest. In present context the only Dirac particles with non-zero bare mass 
are the top quarks and the charginos, whose right- and left-handed plasma masses are 
equal. The bare mass of the charged gauginos W ± is the SU{2) L soft supersymmetry 
breaking gaugino mass M 2 and the bare mass n of charged Higgsinos if ± emerges from 
the bilinear term 11H1H2 in the superpotential. 

We assume that the only particles light enough to be in equilibrium with the 
thermal bath at temperature T are the SM particles, the right-handed stop t R , the 
charginos W ± and H ± , the neutralinos B, W3 and H® 2 , and the neutral and charged 
scalar fields in the two Higgs doublets. This amounts to assuming that all the other 
super symmetric particles, e.g. left-handed stops and gluinos, are much heavier than T. 
This choice is motivated by considerations about the strength of the phase transition 
[0, H |8[ and the mechanism for baryon number generation || f|, |5| in the MSSM. 
For instance, a strongly first order electroweak phase transition can be achieved in 
the presence of a top squark not much heavier than the top quark [7], [J. In order to 
naturally suppress its contribution to the parameter Ap and hence preserve a good 
agreement with the precision measurements at LEP, it should be mainly right-handed. 
This can be achieved if the soft supersymmetry breaking mass nig of is much larger 
than Mz so that ti is decoupled from the thermal bath. It is important to keep 
in mind that this mass hierarchy between the left- and right-handed stops may be 
relaxed if the strength of the phase transition is enhanced by light boson degrees of 
freedom other than the right-handed stops themselves. This happens, for instance, if 
the MSSM content is increased by adding a standard model gauge singlet. 

The relevant hight T plasma masses are given in the Table 1. We have separated 
the SM contribution m§M from the MSSM contribution h-mssm? so that the plasma 
mass is actually the sum of the two. We also denote by g\ = 0.247, g 2 = 0.640 
and #3 = 1.243 respectively the values of U(l)y, SU{2)l and SU(3) C couplings at 
T ~ Mw, and we take the (top-quark) Yukawa coupling to be h t = 1. Other couplings 
and interactions are not included because of their relative smallness. The relevant 
fermionic one-loop diagrams are those, where the right-handed stop couples to one 
of the pairs H^l, H 2 b L or Bt R . The possible bosonic one-loop diagrams consist of 
those including t R and U(l) gauge boson B or t R and gluon. 

Let us first consider the fermionic contribution to the one-loop imaginary part of 
the right-handed stop self-energy. The right-handed stop can absorb a right-handed 
Higgsino hole (left-handed quark hole) to produce a left-handed quark particle (right- 
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Table 1: Plasma masses of light particles contributing to stop decay rate. 



handed Higgsino particle), or it can decay into a pair of particles or a pair of holes. 
The Higgsino may be either a neutralino or a chargino. In the high temperature 
approximation, where the bare masses of the particles running in the loop can always be 
neglected but hard thermal corrections to the propagators are included, the absorption 
contribution to Im E# is given by [1T2J 
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Here (R,h) and (L,p) refer to right-handed holes and left-handed particles, respec- 
tively, and u>L, p (k) and ujr^ are solutions to the fermion dispersion relations, given 
by 
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and H. The coupling factor e 2 f = h 2 for the loops involving and bi (at high T their 



where the upper and lower signs refer to particles and holes, respectively, and m is 
the appropriate plasma mass. The solutions are depicted in Fig. p] for the case of t 
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Figure 1: Fermionic dispersion relations: (a) t-hole, f/2-particle; (b) t-particle, H 2 - 
hole. 



Taking t R to be at rest, the momentum k is determined by the energy conservation 
condition 

M iR + uJ R:h = uj l , p (M ia + ui L>h = uj R:P ) . (9) 



At high T the direct decay channel contributes an imaginary part [PJ 
Ae)k\ul p -k*){u% p -k^ - 
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l-n + ((3u R , p )-n + {(3uj L , v ) (10) 



where k is determined by 



Ml „=UR, a + U L , a (11) 



with a — p,h. At high temperatures the plasma mass of t R is M^ R {T) = 1.020T, 
whereas M^(T) = 0.329T, M$(T) = 0.116T, M H (T) = M tL (T) = 0.649T, and 
Mt R (T) = 0.625T so that decay is clearly kinematically possible to both Bt R and 
H 2 tL or H^l- The full fermionic high T contribution to Im is thus given by the 
sum 2 Im E abs (H 2 t L ) + Im E abs (Bt R ) + 2 Im E deca2/ (^ 2 ^) + Im E decay (Bt R ), and it can 
only be computed numerically by solving k from Eqs. (fl, [IT]) for a fixed Mj (T), using 
the fermionic dispersion relations for holes and particles as given by Eq. (|7|). (The 
factor 2 comes about because the equality of H 2 ti, and H 2 bi contributions.) 

The high T fermionic contributions to T R are depicted in Fig. |2|, plotted against 
Mf R (T), which is related to the zero temperature stop mass by (T) = M| (0) + 
(f #3 + f^f + Wi9i) T 2 - Note that for the fermionic contributions the high T approxi- 
mation does not concern t R ; the results apply equally to the case where t R is no longer 
relativistic (T <C as well as close to the phase transition, where cancellation 

between the zero temperature and high temperature mass terms is possible so that 
T ^> Mi R . Thus for a fixed stop mass, from the figures one can read the various 
contributions to Y R at a fixed temperature. 
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Figure 2: Fermionic contributions to Tr at high T: (a) B-t R -contribution; (b) H-ti 
-contribution. The curves refer to absorption and decay channels as explained in the 
text. Both the decay rate and the mass are given in the units of temperature in a 
logarithmic scale. 

The one-loop, high T gauge boson contribution to Im E# can be cast in the form 



12 



Im X 



e 2 u> g M? R k* 



n {f3u)t) - n (f3u g 



vr n(M 2 - U)ui 

where H(uj g , k) is the longitudinal gauge boson self energy, given by 
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uj g is solved implicitly through u 2 = k 2 + U, and uf = k 2 + M~ r . The coupling factor e 2 
are 4^1/3 and 4g 2 /9 for gluons and B, respectively. The energy conservation condition 
in this case is 

(14) 
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9 ■ 



Physically Eq. (p~2[ ) corresponds to absorption; for kinematic reasons decay is not 



possible for gauge loops. Here we have implicitly assumed that T ^> [Mf (0)| so 
that Mf R (T) < 1.02T. (If we assume M 2 ^(0) to be negative, it is possible that 
Mr (T) T.) This is necessary because of the implementation of the high T ap- 
proximation in the loop, which now also concerns tn as it is circulating in the loop. 
Hence the gauge boson contribution Eq. ( |T2"D is not correct for non-relativistic tn, but 
it is nevertheless valid even if there is a cancellation between the stop bare mass and 
plasma mass terms. 

The gauge boson contributions to T R at high T, plotted against (T), are shown 
in Fig. |3|. One sees that gauge bosons do not contribute to T R when i R is in full 
equilibrium, but may be important if Mf (T) < T. In Fig. ^ we also show the total 
high T thermalization rate T R . At high temperature, far from the the criticial point, 
one finds that at one loop ~ 8 x 10 _4 T. 
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Figure 3: (a) The gauge boson contributions to Tr at high T; (b) the total Tr. Both 
the decay rate and the mass are given in the units of temperature in a logarithmic 
scale. 

The total one-loop result presented in Fig. |3| has a range in which Tr is vanishingly 
small. This reflects only the kinematics of absorption and decay and will disappear 
at higher loops. For instance, at two loops the cuts through various diagrams give 
rise in addition to absorption and decay (involving more than two particles in the 
initial or final state) to processes which correspond to ordinary scattering. An exam- 
ple of the two-loop contributions are the so-called sunrise- diagrams appearing both 
from gluon coupling and quartic Yukawa coupling (the electroweak gauge boson can 
be neglected as discussed earlier). Neglecting the mass differences between different 
loop particles, the imaginary parts can be estimated as in |l2f . We obtain T^ht — 
h 4 t (T 2 /1287rM iR (T)) = 0.0024T 2 /(M fii (T)) and T^ gluon ~ 7^T 2 /(5767rM ffl (T)) = 
0.0092(T 2 /Mr (T)). There appear also t-channel processes from gluon interactions. 
These may be also important whenever the right-handed stop is light. By making 

[0 



we 



a rough approximation of the t-channel gluon processes in the spirit of ref. 
conclude that the t-channel contribution is comparable to the sunrise-diagrams; we 
estimate that T^scatt — 3 x 10~ 3 T when Mr{T) ~ T. Thus two-loop diagrams are 
expected to dominate the thermalization rate when M^ R (T) < 0.1T, or in the range 
0.5T < M~ tR (T) < T, with Y R ~ l^T 2 /M iR (T). In the range 0.1T < M~ tR (T) < 0.5T 
absorption of gauge bosons dominates the thermalization rate, and it can be as high 
as 0.1T. 

Note that right-handed stop plasma masses with Mf (T) <C T can occur only if 
M~ r (0) is negative. Even though present experimental bounds on the stop masses 
do not exclude such a possibility, one should keep in mind that | Mr (0)| is bounded 
from above from considerations about color breaking at zero temperature 0. If we 
restrict ourselves to electroweak baryogenesis and therefore assume that the critical 
temperature is T c = (9(100) GeV, it is easy to show that Mr (T c ) must be larger than 



about 0.55 T c for masses of the lightest Higgs boson around 80 GeV. (For lighter Higgs 
boson the bound for Mr is even larger.) 

Let us briefly discuss the implications of our findings. As we already mentioned 
in the introduction, the precise knowledge of the thermalization rate of the supersym- 
metric particles is a key ingredient for the computation of the final baryon asymmetry. 
Sizeable decay rates of the particles propagating in the plasma destroy the quantum 
interference out of which the the CP-violating sources are built up and therefore re- 
duce the baryon asymmetry. Small decay rates, on the other side, are relevant when 
the particles reflecting off the advancing bubble wall have comparable masses and res- 
onance effects show up In such a case, the thermalization rates provide the natural 
width of these resonances and as the present calculation demonstrates, in supersym- 
metric theories these depend in a complicated way on the particles involved and their 
plasma masses. 

Even though it is presently believed that the right-handed stops do not play a 
leading role in generating the baryon asymmetry, it is important to emphasize that 
this is only true in the context of the MSSM. There the phase transition is made 
strong by the infrared effects in the right-handed stop sector and the baryon number 
is mainly generated by charginos and neutralinos. Should the source of the strength 
of the phase transition reside somewhere else, i.e. in Standard Model gauge singlets, 
the role of right-handed stops will be comparable to the one played by charginos and 
higgsinos and the knowledge of the thermalization rate Tr is important to obtain a 
precise estimate of the final baryon asymmetry. 

If the stop is non-relativistic with M^ R (T) > T, thermalization is dictated by the 
one-loop thermal decay rate which can be larger than T. Thus in any case the thermal- 
ization of tR is rather fast, as can be seen in Fig. |3|. Since the zero temperature limit 
on (0)| from color breaking || suggests that Mf (T) > 0.55 T, we may conclude 
that during baryogenesis the thermalization rate of a relativistic right-handed stop is 
dominated by two-loop effects (i.e. scattering) with Tr ~ 10~ 3 T. This is so in partic- 
ular because the absorption channels close at Mf (T) ~ 0.51 T before decay channels 
open at M$ (T) ~ 0.66 T. This means that the processes of quantum intereference, 
necessary build up the axial stop number, may be damped by the incoherent nature of 
the plasma if tR is non-relativistic at baryogenesis. Even in the favorite case in which 
the left-handed and right-handed stops have comparable masses, the resonance effects 
will be washed out by the large right-handed stop thermalization rate. 

We hope to present our results about the thermalization rate of charginos and 
neutralinos soon flTl |. 
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